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Certain geometric properties of submanifolds of configu- 
ration space are numerically investigated for classical lattice 
cp 4 models in one and two dimensions. Peculiar behaviors of 
the computed geometric quantities are found only in the two- 
dimensional case, when a phase transition is present. The 
observed phenomenology strongly supports, though in an in- 
direct way, a recently proposed topological conjecture about 
a topology change of the configuration space submanifolds as 
counterpart of a phase transition. 
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In Statistical Mechanics phase transitions are asso- 
ciated with the appearence of the so-called Yang-Lee 
(real)zeros |jj of the grand-partition function entailing 
singular temperature-dependence of thermodynamic ob- 
servables. However in the Yang-Lee theory the necessary 
conditions for the existence of real zeros remain unspec- 
ified. It has been recently suggested that the ther- 
modynamic singularities might have their deep origin in 
some major topological change in configuration space, i.e. 
in a non-trivial structure of the support of the equilibrium 
statistical measure. 

This topological conjecture [^j has been put forward 
heuristically within the framework of a numerical in- 
vestigation of the Hamiltonian dynamical counterpart of 
phase transitions. An interesting outcome of such investi- 
gations was the clear evidence of a peculiar temperature- 
behavior of the largest Lyapunov exponent at the phase 
transition point. This was observed in lattice scalar and 
vector ip 4 models |||] , in two and three dimensional XY 
models @ , in the 0-transition of homopolymeric chains 
Pi, and - analytically - in the mean- field XY model §. 
Moreover, in the light of a Riemannian geometrization 
of Hamiltonian dynamics - where Lyapunov exponents 
are related with average curvature properties of subman- 
ifolds of configuration space [p|~|lO|] - also the tempera- 
ture dependence of abstract geometric observables has 
been investigated. Among them the averages of curva- 
ture fluctuations (that enter the analytic formula for the 
largest Lyapunov exponent fiofl ) exhibit a cusp-like pat- 



tern. The peak coincides with the phase transition point. 
Qualitatively similar peaks of curvature fluctuations have 
been reproduced in abstract geometric models (families 
of surfaces of R 3 where the variation of a parameter leads 
to a change of the topological genus) , whence the heuris- 
tic argument about the topological meaning of the peaks 
of configuration-space-curvature fluctuations at a phase 
transition. These fluctuations have been obtained as time 
averages, computed along the dynamical trajectories of 
the Hamiltonian systems under investigation. Now, time 
averages of geometric observables are usually found in 
excellent agreement with ensemble averages P,|3|, ^0| , p| 
therefore one could argue that the mentioned singular- 
like patterns of the averages of geometric observables are 
simply the precursors of truly singular patterns due to 
the tendency of the measures of all the statistical ensem- 
bles to become singular in the limit N — > oo when a phase 
transition is present. In other words, geometric observ- 
ables, like any other "honest" observable, already at fi- 
nite N would feel the eventually singular character of the 
statistical measures, and if this was the true explanation 
we could not attribute the cusp- like patterns of curvature 
fluctuations to special geometric features of configuration 
space. Hence the motivation of the present paper. We 
aim at elucidating this important point by working out 
purely geometric informations about the submanifolds (to 
be specified below) of configuration space, independently 
of the statistical measures. In the present paper a new 
step forward is done in the direction of supporting the 
topological conjecture mentioned above. 

Our geometrical framework is the configuration space 
M of systems whose degrees of freedom are uncostrained 
and are real numbers, thus M = K . 

Our statistical framework is the canonical ensemble 
whose volume in M is given by the configurational par- 
tition function Zq — f YliLi dqi exp[—V (q)] , where q = 
(qi, . . . , 3jv) S M. . Finally, our topological framework is 
elementary Morse theory JlJ|l3|]. Let us here recall its 
basic idea with the help of Fig. 1. The "U-shaped" cylin- 
der of Fig.l is the ambient manifold M where a function 
V : M i— ► M, smooth and bounded below, is defined to be 
the height of any point of M with respect to the floor- 
plane. For any given value u of the function V two kinds 
of submanifolds are determined: the level sets S u of all 
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the points x S M for which V(x) = u, and M u , the part 
of M below the level u, i.e. the set of all points x 6 M 
such that V(x) < u. The remarkable fact of Morse the- 
ory is that from the knowledge of all the critical points 
of V, i.e. those points where W = 0, and of their Morse 
indexes, i.e. the number of negative eigenvalues of the 
Hessian of V, one can infer the topological structure of 
the manifolds M u , provided that the critical points are 
nondegenerate, i.e. with nonvanishing eigenvalues of the 
Hessian of V. Two such points are marked in Fig.l, at 
the bottom of M and at some intermediate level u c for 
which S Uc is an eight-shaped curve. It is evident from 
Fig.l that the manifolds M u<Uc are not diffeomorphic 
to the manifolds M u>Uc : the formers are homeomorphic 
to a disk and the latters are homeomorphic to a cylin- 
der. The same happens (in general) to the boundaries S u 
that here are circles for u < u c and become the topolog- 
ical union of two circles for u > u c . This simple example 
displays the general fact that passing a critical level of a 
Morse function is in one-to-one relation with a topology 
change. A critical level is a surface E u that contains one 
or more critical points. 

Let us now consider the configuration space M — R N 
of a physical system and its potential V as the Morse 
function. The interesting things are supposed to oc- 
cur below some large value u so that the corresponding 
(large) subset M C M is compact. Then M and all its 
submanifolds M u are given a Riemannian metric g. On 
all these manifolds (M u ,g) there is a standard invariant 
volume element: dij — y det(g)dq 1 . . . dq N . 

In order to study the topology of the family {(M u , g)} 
we should find, analytically or numerically, all the crit- 
ical points of V, but at large N this is a formidable 
task, therefore we have approached this problem as fol- 
lows. Generalizing a simple geometric example reported 
in we have computed the total degree of second 

variation u 2 K of the Gaussian curvature, i.e. a 2 K = 
(Kq)j: u — (Kg)z where (•) stands for integration over 
the surface £„, as a function of u in the neighborhood 
of a critical point. This is possible in general, indepen- 
dently of the potential V, because any Morse function 
can be parametrized in the neighborhood of a critical 
point, located at xq s M, by means of the so-called 
Morse chart, i.e. a system of local coordinates {xi} such 
that V(x) = V(x ) - + El k+ i4 (k is the 

Morse index). Then standard formulae for Gauss cur- 
vature of hypersurfaces of R N [Q can be used to ex- 
plicitly work out Kg and o~ 2 K . The intersection of an 
hypersphere of unit radius - centered around u = 
(the critical point) - with each S M is used to bound 
the domain of integration. The numerically tabulated 
results are reported in Fig. 2 and show that o 2 K devel- 
ops a sharp, singular peak as the critical surface is ap- 
proached. It seems therefore reasonable to apply this 
geometric probing of the presence of critical points, and 



hence of topology changes, to the study of possible topol- 
ogy changes of the manifolds (M u ,g). In fact the mani- 
folds M u inheritate - by somewhat smoothing them - the 
peculiar geometric properties of the {T> u } u < Uc - due to 
the presence of critical points - because of the relation 
Jm fdv — Jo dv J s /|s v dw/|| VV||, where du> is the in- 
duced measure on £„ and / a generic function jD| . The 
surface S Uc defined by V(x) — u c is a degenerate quadric, 
therefore close to it some of the principal curvatures 0] 
of the surfaces E u ~ Uc tend to diverge. Such a divergence 
is generically detected by any function of the principal 
curvatures and thus, for practical computational reasons, 
instead of Gauss curvature (which is the product of all 
the principal curvatures) we shall consider the total sec- 
ond variation of the scalar curvature 1Z (i.e. the sum of 
all the possible products of two principal curvatures) of 
the manifolds [M u ,g), according to the definition 

a 2 n {u) = [VoliMu)]- 1 [ dr)(K- g u f (1) 
Jm u 

Qu = [Vol{M u )Y 1 [ d V K (2) 

with 1Z = g k: >R l kl j, where R l ki - is the Riemann curvature 
tensor @, and Vol(M u ) = f drj. 

Let us now consider the family of submanifolds 
(M u ,g) associated with the so-called ip 4 model, on a d- 
dimensional lattice Z d , described by the potential func- 
tion 

where (af3) stands for nearest-neighboring sites. We con- 
sider d = 1,2; this system has a discrete Z2-symmetry 
and short-range interactions, therefore in d = 1 there is 
no phase transition whereas in d = 2 there is a symmetry- 
breaking transition (this system has the same universal- 
ity class of the 2d Ising model). 

The potential in Eq. (^) defines the subsets M u of con- 
figuration space, these subsets are given the structure of 
Riemannian manifolds (M ui g) by endowing all of them 
with the same metric tensor g. However, as far as we 
want to learn something about the topology of these 
manifolds, the choice of the metric g is arbitrary. We 
have therefore chosen three different types of metrics, one 
conformally flat and the others non-conformal, according 
to a compromise between simplicity and non-triviality. 
They are: i) g^J = [A — V{q)]5 liv , i.e. a conformal de- 
formation of the euclidean flat metric 5^, A > is an 

(2) (3) 

arbitrary constant larger than u; ii) g^J and g^J are 
non-conformal metrics defined by 

//{*) o 1\ 
(«#?)= 10, fe = 2,3 (4) 
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where I is the N — 2 dimensional identity matrix, </ 2 ) is 
obtained by setting /< 2 ) = ± J2 a ez" it + A > and 5 (3) b y 
setting /( 3 ) = -k J2 aeZ d q^+A, with A > 0, and a labels 
the JV lattice sites of a linear chain (d = 1) or of a square 
lattice (d = 2, N = nx n). Simple algebra Jl6| yields the 
scalar curvature function for each metric: 



n {1) = (N - 1) 



AV 



|W|| 



(A - V) 2 {A - Vf 
||V/Wf 



N 

T 



(5) 



(/(*) - 1) 



2(/W - 1) 



A/ 



(fe) 



fc = 2,3 (6) 



where V and A are euclidean gradient and laplacian re- 
spectively; V and A do not contain the derivatives d/dq a 
with a = 1 (d = 1) or a = (1, 1) (d = 2). 

We constructed an ad-hoc MonteCarlo algorithm to 
sample the geometric measure drj by means of the stan- 
dard "importance sampling" method |l7|j , then we ap- 
plied it to the computation of er^(it), given by Eq.(||), 
for the one- and two- dimensional lattice ip 4 model de- 
fined in (^) with the following choice of the parameters: 
A = 0.6, fi 2 — 2, J = 1. The values of TZ are com- 
puted according to Eqs.(||) and (||). In order to locate the 
phase transition that occurs in the two-dimensional case, 
we have computed (V) vs T by means of both Monte- 
Carlo averaging with the canonical configurational mea- 
sure, and Hamiltonian dynamics (by adding to V a stan- 
dard kinetic energy term). In the latter case the temper- 
ature T is given by the average kinetic energy per degree 
of freedom, and (V) is obtained as time average. Fig. 3 
shows a perfect agreement between time and ensemble 
averages, thus we worked out Fig. 3 by computing 200 
time averages because they converge much faster than 
ensemble ones. The phase transition point is well visible 
at u c = (V) ~ 3.75. 

In Figs. 4 and 5 we synoptically report the patterns of 
a Ti( u ) f° r the one and two dimensional cases obtained 
at different lattice size with g^ 1 ' (Fig. 4), and obtained at 
given lattice size with g( 2 * 3 ) (Fig. 5). Peaks of cr^(u) ap- 
pear at u c - the value of (V) at the phase transition point 
- in the two-dimensional case, whereas only monotonic 
patterns are found in the one-dimensional case, where no 
phase transition is present. 

"Singular" , cuspy patterns of a^(u) (with the meaning 
that such attributes can have for numerical results) are 
found independently of any possible statistical mechani- 
cal effect, and independently of the geometric structure 
given to the family {M u } by the metric tensors chosen. 
Though within the well evident limits of numerical simu- 
lations and of a limited choice of different metrics, our re- 
sults suggest that the "singular" patterns are most likely 
to have their origin at a deeper level than the geomet- 
ric one, i.e. at the topological level. Hence the observed 
phenomenology strongly hints at the occurrence of some 
major change in the topology of the configuration-space- 
submanifolds {M u } in correspondence with a second- 



order phase transition. Finally, the expression " major 
topology change" is to suggest that a change of the co- 
homological type of the M u - or £„ - might well be a 
necessary but not sufficient condition for a phase tran- 
sition, and that in any case some "big" change has to 
occur. 
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FIG. 1. Illustration of the relationship between topology 
and critical points. The U-shaped manifold M is born at Pq. 
The level surfaces E„ and the parts of M below them - M u 
- change topology when u exceeds the height of the critical 
point Pi. 
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FIG. 2. Variance of Gauss curvature vs u close to a critical 
point. a 2 J N is reported because it is dimensionally homoge- 
neous to the scalar curvature. Here N = dim(E u ) = 100, 
and Morse indexes are: k = 1, 15, 33, 48, represented by solid, 
dotted, dashed, long-dashed lines respectively. 

FIG. 3. Average potential energy vs temperature for the 
2d lattice tp 4 model. Lattice size N — 20 x 20. The solid 
line, made out of 200 points, refers to time averages. Full 
circles represent MonteCarlo estimates of canonical ensemble 
averages. The dotted lines locate the phase transition. 

FIG. 4. Variance of the scalar curvature of M u vs u com- 
puted with the metric g^K Full circles correspond to the 
ld-ip 4 model with N = 400. Open circles refer to the 2d-ip 4 
model with N = 20 x 20 lattice sites, and full triangles refer 
to 40 x 40 lattice sites (whose values are rescaled for graphic 
reasons) . 

FIG. 5. Variance of the scalar curvature of M u vs u com- 
puted for the ip 4 model with: metric in lei, N = 400 
(open triangles); metric g^ in 2c?, N = 20 x 20 (full trian- 
gles); metric g^ in lei, iV = 400 (open circles); metric g^ in 
2d, N = 20 x 20 (full circles). 
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Fig. 1 Franzosi et al. 
"Topological aspects of ..." 




Fig. 2 Franzosi et al. 
"Topological aspects of ..." 




Fig. 3 Franzosi et al. 
"Topological aspects of ..." 




Fig. 4 Franzosi et al. 
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